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We discuss the cosmological reconstruction in modified Gauss-Bonnet and F(R) gravities. Two 
alternative representations of the action (with and without auxiliary scalar) are considered. The 
approximate description of deceleration-acceleration transition cosmologies is reconstructed. It is 
shown that cosmological solution containing Big Bang and Big Rip singularities may be recon- 
structed only using the representation with the auxiliary field. The analytical description of the 
deceleration-acceleration transition cosmology in modified Gauss-Bonnet gravity is demonstrated to 
be impossible at sufficiently general conditions. 

PACS numbers: 95.36.+X, 98.80.Cq, 04.50.Kd, ll.10.Kk, ll.25.-w 

I. INTRODUCTION 

The modified gravity approach (for general review, see []]]) became the essential element of the modern cosmology. 
It is quite remarkable that some change of the classical gravitational action may resolve the number of cosmological 
problems, including inflationary paradigm, dark energy and dark matter. It turns out that it is not necessary to 
introduce the extra ingredients (usually, scalar or fluid) as all these phenomena could be understood as gravitational 
manifestations. For instance, the unification of the early-time inflation and late-time acceleration may be achieved in 
F(R) gravity (for first realistic model of that sort, see 01) without the need to introduce the inflaton and (scalar) dark 
energy by hands. Several models of modified gravity may successfully describe dark matter as gravitational effect (for 
a recent review, see [3]). The coincidence problem effectively disappears in the modified gravity approach because 
dark matter and dark energy are caused by the universe expansion governed by specific theory. It is expected that 
modified gravity may be helpful also in high-energy physics (for instance, for hierarchy problem). 

Unfortunately, the realistic modified gravity has usually highly non-linear structure in terms of geometric invariants 
(curvature, Gauss-Bonnet invariant, etc.). As the result, its background evolution is very hard to describe analytically 
unlike to the case of General Relativity where number of viable analytic solutions are available. In turn, with only 
approximate FRW solutions of modified gravity it is extremely difficult to study the cosmological perturbations. At 
best, such cosmological perturbations are studied in further approximation neglecting the higher-derivatives non- 
linearities which is definitely not sufficient. In order to study the background evolution of the alternative gravities, 
so-called reconstruction method has been developed (for the introduction, see [![)■ Within the reconstruction method, 
given FRW cosmology may be used to reconstruct the modified gravity where such cosmology is the solution of the 
equations of motion. 

In the present paper we develop the reconstruction method for modified Gauss-Bonnet gravity [5| . It is demonstrated 
how to reconstruct the theory which admits the deceleration-acceleration transition (the transition to ACDM epoch) . 
Such background evolution turns out to be the very complicated and approximate one. For quite general class of 
_F(G)-functions we show that there is no analytical description of deceleration-acceleration transition. It turns out 
that it is very difficult (if possible at all) to construct such a model which admits such transition analytically (the 
co-existence of matter dominance and accelerating solutions [6]). The comparison with F(R) — R + f(R) theory is 
done. The alternative presentation for F(R) and F(G) modified gravity using the auxiliary scalar is considered. It 
is shown that reconstruction using such representation leads to wider class of cosmological solutions, including the 
deceleration-acceleration transition ones. It is demonstrated that cosmological solution containing the Big Bang as 
well as Big Rip singularity may be reconstructed from F(R) gravity. 
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II. ANALYTICAL APPROACH TO DECELERATION-ACCELERATION TRANSITION IN 

/(G)-GRAVITY 



Let us study modified gravity with the following action 



S F (G) = J d^Xy/^g (^J^ + F(G) + Anattcr^ 



(1) 



Here G is the Gauss-Bonnet invariant G — R 2 — AR^R^ + R l _ u ^ p(T R >i ^ f " 7 and £ ma tter is the Lagrangian of matter. It 
is convenient to put 2k 2 = 1 in this section. We will discuss only the FRW background: 

g^ v = diag(-n 2 ,a 2 ,a 2 ,a 2 ) . (2) 
The variation of the action (Q]) with respect to lapse-function n gives the modified Friedman equation: 

6H 2 + F- F'G + 24F"H 3 G = p, (3) 
and the variation with respect to scale factor a gives more complicated equation: 

2C 

4H + 6H 2 + F- F'G + 8F"'H 2 G 2 + —^F"G + 8H 2 F"G = -p . (4) 

Here p and p is the matter energy- density and pressure, respectively, which arises from £ ma ttcr- We also note that prime 
' denotes partial differentiation of function F with respect to its argument. Using analogy with the FRW equations in 
the Einstein gravity one may define p G = —F + F'G—2AF"H a G and pa = F-F'G + 8F"'H 2 G 2 + ^F"G + 8H 2 F"G, 
so the equations ^ and Q take the following form: 

6H 2 = Ptot , (5) 



AH + 6H 2 = - Ptot , (6) 

where p to t = P + Pg an d Ptot = P+Pg- Note that different cosmological solutions for above theory have been discussed 
in refs.Q. 

The barotropic equation of matter state p = wp is considered below. To close our system one adds conservation 
energy law which takes the following form: 

p + 'SH(p + p)=0. (7) 

Note also that equation ([3]) is just the first integral of the system (@J-([7]). Now let us consider the possibility of 
occurrence of late-time universe acceleration due to function F(G) analytically. In other words, one searches some 
function F(G), which plays non-essential role during the standard dust stage(w = 0), but gives leading contribution at 
late times. Our purpose is the analytical description of such deceleration-acceleration transition. We discuss functions 
with F(0) = 0, because otherwise we will have some analogue of cosmological constant. (Of course, permitting the 
effective cosmological constant may qualitatively change the results obtained below) . The effective equation of state 

parameter may be easily found by using the expressions (JS])-®: w c s = — 1 — j§z- 

Now let us suppose that there exists some function F(G) which leads to late-time acceleration and has the following 
properties: 

F(0) = 0, F'(0) = 0, F"(Q) = 0, F'"(0) = . (8) 

From another side it is known that at the deceleration-acceleration transition point G — 24^# = 0. This point is 
reached when w to t = —1/3. By calculating the values of p G and pq at the transition point, we find that they vanish. 
This means that there is no any effective matter besides the usual matter p at this moment, so ptot = P and ptot = P, 
hence i«tot must be equal to some value of w which is bigger than —1/3. This logical contradiction proves that any 
function satisfying the conditions ([5]) cannot reach the deceleration-acceleration transition point. Note also that the 
condition ^'(0) = may be removed from (|8]) because its contribution to pa and pa contains G as a factor. This 
result is complimentary to the one of ref. where it has been shown that some class of F(G)-theories which allow an 
exact power-law solution can not explain transition from deceleration to acceleration. Actually, the following general 
form of function which allows exact power-law decelerating solution (see (15) in [6]) is: F{G) — AG 5 + BG k where 
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k < | and this function does not satisfy to our condition ((SJ). From another side, one may easily find functions which 
satisfy the conditions ([5} but do not allow exact power-law solution. For example, 

F = J> G », r-^r^. (9) 

N=4 



The situation is the following: there is some function F, which does not allow exact power-law decelerating solution, 
but allows it approximately with very good accuracy. This solution may be unstable and leading to acceleration. The 
examples of such approximate deceleration-acceleration transition will be discussed below. 

III. COMPARISON WITH F(i?)-GRAVITY 

It is interesting to compare results of the previous section with .F(i?)-gravity. Its action has the following form: 



S m) = J a^x^-g ( JL [R + f(R)} + £ mattOT ) 



(10) 



FRW equations of motion are (here again 2k 2 = 1): 

6H 2 + f - f'R + 6H 2 f + 6Hf"R = p, (11) 

and 

AH + 6H 2 + / - f'R + 6H 2 f + 2f"'R 2 + AH f'R + 2f"R = -p . (12) 

The latter equation is a consequence of (|TT1) and ([7]). It is well known that R = identically for the regime a ~ i 1 / 2 , 
which corresponds to w c tf = |. So using developed analysis one may try to investigate the possibility to reach a ~ t 1 / 2 
regime. We will study only functions satisfying the conditions 

/(0) = 0, /' (0) = 0, /" (0) = 0, f"'(Q) = . (13) 

Let us consider the theory with fixed EoS w m matter besides the relativistic matter (w m — ^). Following previous 

section we define p f(R) = -f + f'R - 6H 2 f - QHf'R and p f{R) =f- f'R + 6H 2 f + 2f"'R 2 + AHf'R + 2f"R to 
rewrite equations (fTTj) -(|12 p in the canonical form In this case we have the logical contradiction: from the one 

side it must be Wtot = ^ on a ~ t 1 / 2 regime, but from another side we have wtot = w m ^ ^ b ecause there is no any 
contribution to /9 to t and p t ot from /(i?)-terms at this regime (pf(R) = 0, Pj(r) — due to (|13|l ). 

So we have the following result. The analytical description of transition to regime a ~ t 1 / 2 in the universe with 
any perfect fluid except w m — ^ in /(i?)-gravity with p3p is very hard to realize (compare with [8| where similar 
conclusion is made). Of course, other classes of functions f(R) or account of the effective cosmological constant may 
change this conclusion. 

Note also that there is no any problem with deceleration- acceleration transition in f(R) gravity. A number of 
such theories admitting the transition is well known. For example, most general function which leads from matter 
dominated era to the ACDM cosmology was constructed in [j| [Tfl] by using reconstruction method. Below we try to 
use this method to solve the problem described in previous section. 



IV. RECONSTRUCTION AND THE DECELERATION-ACCELERATION TRANSITION IN 

F(G)-GRAVITY 

Let us investigate possibility to find the theories ([1} which allow transition from deceleration to acceleration phase by 
using reconstruction method. This method developed in ref.[9( may be easily adopted to our Gauss-Bonnet modified 
gravity ([T]). We start from the equation First of all we would like to use a new variable N instead of the 

cosmological time t, defined by N = In Here ao is the value of the scale factor a(t) in ([2]) at a fixed time. This 

variable is related with the redshift z by e~ N = 1 + z. Since = H-^ and = R 2 -^i + ^^7737^1 one can rewrite 
© as 



6H 2 + F(G) - 24H 3 {H' + H)F'(G) + 24 2 F" (G)H 6 (HH" + 3H' 2 + AHH') = p . 



(14) 
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Here W = dH/dN and H" = d 2 H/dN 2 , but F' = dF/dG like above. Here we have used also G = 2AH 2 {H + H 2 ) = 
24H 3 (H' + H). If the matter energy density p is given by a sum of the fluid densities with constant EoS parameter 
Wi, we find 



-3(1+^.,) _ \p 



-3(l+u)i) -3(l+ Wi )N 



P = Y2pioa " v ~ ■ —/ = y ^pioOp ' ■ "e 

Here pio is a constant. Let the Hubble rate is given in terms of TV via some function k(N) as 

H = k(N) =k(- ln(l + z)). 

Note now that the expression G = 2Ak(N) 3 k' (N) + 24/c(A r ) 4 may be solved with respect to N as N 
by using (fl~5j) and (fl6|) . one can rewrite (|14| as 



(15) 
(16) 

N(G). Then 



6 (Jfe (TV (G))) 2 + F(G) - 24 (k (N (G))) 3 F'(G) \k' (N (G)) + k {N (G))] 
+24 2 F"(G) (jfe (N (G))) 6 [ft (JV (G)) k" [N (G)) + 3 (k' {N (G))) 2 + 4fc [N (G)) fc' (iV (G)) 

This equation is differential equation for F(G) and may be simplified by introducing h(N) = (k (N)) 2 = H 2 : 

dF{G) 



(17) 



6h(N(G)) + F(G) - 12- 
<2d 2 F(G) 



+24 z 



rfG 2 



-fc(JV(G0)- 



{ h(N(G))h'(N(G))+2(h(N(G))y 
\«'(N(G )) + 2k'(N(G)) + ^^ 



(18) 



Note that the Gauss-Bonnet invariant is given by G = 24h(N) 2 + 12h(N)h'(N). Hence, when we find F(G) satisfying 
the differential equation (|18p. such F(G) theory admits the solution ([TBI and therefore such gravity realizes above 
cosmological solution. This is essentially the cosmological reconstruction. 

Now let us discuss the simplest example which is related with previous discussion and which reproduces the ACDM- 
era. In the Einstein gravity the FRW equation for the ACDM cosmology is given by 



6H 2 = 6#o + Poa 



6H 2 



p Q a 3 e m . 



(19) 



Here Hq and po are constants. This equation reproduces the universe with dust matter which enters to ACDM-era 
at late time (for sufficiently small Hq). Therefore, it reaches the point d — at some moment. So we have 

1 



h(N) 



Hi 



Substituting this relation into expression for G we find: 

G(N) = 24ij4 + 2H 2 p a- 3 e- 3N 



± fl 2 -6 -6JV 



which may be solved to find N(G). It is convenient to introduce x 

1 



2Am 



2Hlx 



Po% 
2 =G. 



-3N 



, so finally we have 



(20) 



(21) 



(22) 



It is interesting to note that x > at any moment. Moreover, one may easily calculate the transition point from 
deceleration to acceleration which corresponds to G = 0: x(G = 0) = 12Hq. The solution of (p?2")) is: 



Xl ,2 = 3H 2 ± 



9 2 H$ 



3G. 



(23) 



Note that sign "— " must be excluded because it corresponds to non-physical negative values of x for negative G, which 
corresponds to accelerated regimes. Now we can see that only values G < 27Hq are resolved. Now it is necessary to 
use the function (|20| in order to find the theory F(G) as the solution of the differential equation (fT8|). It turns out that 
this differential equation is extremely complicated and the corresponding solution may be found only numerically for 
different asymptotics (near to transition point). In principle, it is easier to construct such solutions in the alternative 
representation of F(G) theory with auxiliary scalar. The corresponding examples are found in third and fourth papers 
from ref.Q. That is why we will no go further to technical details of the solution of eq. fTg)) . Hence, in principle it is 
possible to construct F(G) which allows the transition from deceleration to acceleration era (of course if we can solve 
the corresponding differential equation). 
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V. ALTERNATIVE REPRESENTATIONS OF F(i?)-GRAVITY AND F(G)-GRAVITY AND THE 

RECONSTRUCTION 



Let us discuss the alternative representation for F(i?)-gravity and F(G)-gravity with the actions given by (fT0|) and 
([T]). respectively. 

In addition to the problems mentioned in the previous sections, there appear other problems in F(R)- and F(G)- 
gravities. First problem is easy to understand in terms of F(R) gravity. The Einstein gravity coupled with perfect 
fluid with constant equation of state (EoS) parameter w can be reproduced by the following F(R) theory 



um 9w + 7± V45w 2 + 126w + 53 2(m-2) 

-F(i?) (x R , m= —. TT or w = -1 - — w „ . (24) 

v ; ' 6(w + l) 3(m-l)(2m-l) v 7 

We may investigate the modified gravity with Big Bang singularity with u>bb > and Big Rip singularity [Tl| with 
w = wbr < 0. In both of the Big Bang singularity and Big Rip singularity, the scalar curvature R diverges. This 
shows that if we construct a model describing both of Big Bang singularity and Big Rip singularity, the corresponding 
F(R) must be double valued function of R. 

Similarly for i 7 '(G)-gravity, if we try to construct a realistic model, where there is a transition from decelerating 
phase to the accelerating phase, F(G) may become a double valued function or it may become purely imaginary 
function. 

In the following, we consider how the above problem could be solved. At least locally we can rewrite the actions 
pop and (JTJ) by introducing the auxiliary scalar field 4> as follows 

Z>F(R) = dx yj-g I — h£matter I , (25) 



and 



R 



S F (G) = / dx 4 V=9 7, o - + /W>) G + ^tter ■ ( 26 ) 



F(G) = I ax V~9\^2 

We should note, however, the actions (f2~51 and (|26p express more wide class of theories than the actions (flU)) and 
(HJ) (for related discussion, see also [12j). For example, we may consider the following model corresponding to F(R) 
gravity: 

P(<t>) = \4? + P4>\ Q(<l>) =i<t>- (27) 

Here (3 and 7 are constants. (The following arguments do apply even for F(G) gravity.) Then by the variation of <fi, 
one finds 

= <P 2 + 2/30 + jR , (28) 
0=-/3± V/3 2 -lR, (29) 



which can be solved with respect <f> as 



which gives 

S F{R) = J dx^^—g + £ matter ) , F ± (R) = + ^ (-p ± ^^R 



(30) 



The action (|30|) is double-valued function and furthermore the value of R is restricted to be jR < 1 in order to have 
the real Sp(R)- Hence, the action (fT0|) describes the theory corresponding to one of the branches of double- valued 
function and R is restricted to be jR < (3 2 . We should note, however, that we need not to start from the action 
(flTJ| but from the action (|25|) . The action (|25jl may describe the scalar field theory with potential — ^jjr and the 

Brans-Dicke non-minimal coupling but without the kinetic term. If we start with the action (|23|) instead of (flU)) 
from the very beginning, even if we consider the model (|2T|) , we may obtain the theory with transition between F + (R) 
and F-(R). Note that, in the model corresponding to (|23|) with (|2"T|) , the value of R can be, in general, in the region 
7-R < /3 2 , which is forbidden for the action (flU)) . 







Let us clarify it in more detail. For simplicity, we neglect the contribution from matter by omitting £ ma tter- First 
we consider the following model corresponding to (|25p: 



P(0) = e^)/ 2 p(0) , g(ct>) 



-10 In 



C 



7+1 



P{4>) = P4 



3 + 



/3± 



i ± Vi + 1007(7 + 1) 



Q(0) 



-G 



P(<£) - 6 



dg{4>) dP(</>) 



Here to, C, and p± are constants. Now P(</>) and Q(</>) are smooth functions of as long as 

< <p < t s = toC-^ + V . 

The exact solution of the FRW equation is 



H(t) 



Hi)' 7 ' 1 + (7 + l)g(^) 7 



(31) 



(32) 



(33) 



When t — > 0, i.e., i <C i s , -ff(i) behaves as 

which corresponds to the Big Bang singularity at t — 0. On the other hand, when t — > t s , we find 

Hit) 10 



t. - < 



(34) 



(35) 



which corresponds to the Big Rip singularity. Then in the form ([23)1 of the action, one can obtain the cosmological 
model describing both of the Big Bang and Big Rip singularities. In this alternative presentation with auxiliary scalar, 
it is easy also to construct the deceleration-acceleration transition solutions. Such reconstruction has been presented 
already for F(R) and F(G) theories in refs.JJ, [l0| and third and fourth papers from ref.[7|. That is why we do not 
give the details of such cosmologically-viable theories here. 

We should also note that the modified gravity which exhibits the transition from deceleration epoch to acceleration 
epoch can be obtained without introducing the auxiliary field (j) (compare with reconstruction in refs.(9l. [ToL fl3l|). For 
example, we consider the following form of F(R): 

F(x) = AF(a, /3, 7; x) + Bx l ~^F[a - 7 + 1, /3 - 7 + 1, 2 - 7; x) . (36) 

Here A and B are constants, F(a,/3,j;x) is Gauss' hypergeometric function, x is defined by x = — 3, and 
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-2' a 



-6' a/? 



(37) 



The action has an exact solution which reproduces, without real matter, the ACDM era whose FRW equation is given 
by 



K K 

Next we consider the following F(G) gravity model corresponding to the action 

2 „ / , X / , x 9l 



(38) 



■>o 



651 



>o 



f(4>) 



o3i f *° dx eX x 91 



4k 2 



(1 + gix) 



4%k 



Hi 



(39) 
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Here g\ and 4>o and positive constants and $W(a,x) is given by the incomplete gamma function: 

W(a,x) = f dye'Vy"- 1 . (40) 



Note that the functions V{<f>) and /((/>) are smooth functions as long as <j> > 0. An exact solution of the model is given 
by 

00 t 

When t is small H(t) describes the Big Bang singularity where the expansion of the universe is decelerating if g\ < 1. 
On the other hand, when t is large H goes to a constant: H —> 4^, which corresponds to the de Sitter universe and the 
universe is expanding with the acceleration.. Hence, starting from the theory with the action (|26[) . one can explicitly 
construct a model which admits the approximate transition from decelerating phase to the accelerating phase, 



VI. DISCUSSION 



In summary, we discussed the cosmological reconstruction method for modified Gauss-Bonnet and F(R) gravities. 
Two alternative representations for the action is used: with and without the auxiliary scalar field. It turns out that 
the cosmological solutions in the representation with the auxiliary scalar follow from the wider class of theories. 
Moreover, it is easier to reconstruct modified gravity in such representation. For instance, the cosmological solution 
which contains the Big Bang and Big Rip singularities may be reconstructed in such formulation with the auxiliary 
scalar but not in the original formulation. Special attention is paid to the cosmologies admitting the deceleration- 
acceleration transitions. It is shown that such cosmological solutions may be reconstructed in both representations of 
modified gravity but only approximately. The analytical deceleration-acceleration transition cosmology in modified 
Gauss-Bonnet gravity satisfying to some reasonable conditions is shown to be impossible. It is extremely hard (if 
possible at all) to find such analytical solutions in modified Gauss-Bonnet gravity. 

The detailed understanding of the background evolution of modified gravity is the necessary step in the development 
of the cosmological perturbations. Hence, even the approximate background evolution realized via the reconstruction 
method may serve for this purpose in order to select the most realistic theories confronting them with the observational 
data. 
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